plu decomposition method